Oscillations of magnetization and conductivity in anisotropic 
Fulde-Ferrell-Larkin-Ovchinnikov superconductors 
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We derive the fluctuational magnetization and the paraconductivity of Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) superconductors in their normal state. The FFLO superconducting fluctu- 
ations induce oscillations of the magnetization between diamagnetism and unusual paramagnetism 
which originates from the competition between paramagnetic and orbital effects. We also predict a 
strong anisotropy of the paraconductivity when the FFLO transition is approached in contrast with 
the case of a uniform BCS state. Finally building a Ginzburg-Levanyuk argument, we demonstrate 
that these fluctuation effects can be safely treated within the Gaussian approximation since the crit- 
ical fluctuations are proeminent only within an experimentally inaccessible temperature interval. 
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I. INTRODUCTION. 

Forty years ago, Fulde and Ferrell^, and Larkin and 
Ovchinnikov^ predicted that the paramagnetism of the 
electron gas might induce a novel superconducting state 
wherein the order parameter is modulated in real-space. 
In their original proposal, these authors considered a sin- 
glet s-wave superconductor perturbed by the Zeeman ef- 
fect (paramagnetic effect), and neglected completely the 
orbital coupling and the disorder. For most type-II su- 
perconductors, the superconductivity is destroyed by the 
orbital pair-breaking effect which leads to a more famil- 
iar inhomogeneous superconducting state: the Abrikosov 
vortex lattice. In order to observe the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state, the paramagnetic ef- 
fect must break Cooper pairs more efficiently than the 
orbital one. Such a situation may be realized in tridi- 
mensional (3D) superconductors with large internal ex- 
change fields, like the rare-earth magnetic superconduc- 
tor ErRh4B4, see"^ for a review. Another possibility cor- 
responds to a quasi two-dimensional (2D) layered super- 
conductor wherein the weakness of the interplane hop- 
ping suppresses the orbital effect for in-plane magnetic 
field. Being the ratio of the critical fields, H°''''{T = 0) 
in the pure orbital limit, and Hp(T = 0) in the pure para- 
magnetic limit, the Maki parameter am — \/2H°^'''{T = 
0)/Hp{T — 0) quantifies the relative strength of those 
pair breaking mechanisms. Besides demanding a large 
Maki parameter (a,„ > 1.8), the occurance of the FFLO 
state also requires very clean samples since it is far less 
robust against disorder than the usual vortex lattice, see"* 
and^ for recent reviews. 

Recently, there have been mounting evidence that the 
heavy fermion superconductor CeCoIn5 under magnetic 
field might fuUfil those stringents conditions'"^, although 
the magnetism of this system is still under debate. A 
superconducting phase have been reported at large mag- 
netic field and low temperature which is distinct from the 
uniform superconducting phase realized at lower fields. 
The characteristics of this phase depends upon the ori- 
entation of the field relatively to the basal plane of the 
tetragonal CeCoIns lattice^ . In the field- induced organic 



superconductor^ A-(BETS)2FeCl4, and in the layered or- 
ganic superconductorio k-(BEDT-TTF)2Cu(NCS)2 the 
FFLO state have been reported when a strong magnetic 
field (20 T for latter one) is applied along the supercon- 
ducting planes. 

However in practice, the identification of the FFLO 
state is hindered by the interplay between orbital and 
paramagnetic effects. The first available experimental 
clue is the shape of the transition line Hc{T) separating 
the normal state from the inhomogeneous superconduct- 
ing state. A lot of theoretical works have been devoted 
to the description of this Hc{T) line. For moderate Maki 
parameters, am < 9, the structure of the FFLO modula- 
tion involves a zero Landau level (index n — 0) function 
(Gaussian with no additional modulation) For higher 
Maki parameter, am > 9, the Cooper pair wave function 
of a 3D superconductor consists in a cascade of more ex- 
otic solutions, the so-called multi-quanta states, which 
are described by a higher (index n > 0) Landau level^^. 
Such values of Maki parameters are rather high for 3D 
compounds (for instance CeCoIns has am = 4.6 — 5) but 
they can be achieved in layered quasi-2D superconductors 
(or superconducting thin films) under in-plane magnetic 
fields^'^. All these studies were performed so far in the 
framework of isotropic models, namely for the idealis- 
tic case of a spherical Fermi surface in the normal state. 
Moreover it has been shown that an elliptic Fermi surface 
leads to the same phenomenology at cost of introducing 
an angle-dependent Maki parameter^'*. 

In real compounds, the crystal lattice (or the pairing 
symmetry) induces a non trivial anisotropy which mat- 
ters a lot for the modulated state^^'^' since it essentially 
lifts the degeneracy between various orientations of the 
FFLO modulation. Recently, the interplay of paramag- 
netic and orbital effects was reconsidered in the presence 
of such a non trivial anisotropy, namely for a Fermi sur- 
face which slightly differs from the spherical or elliptical 
shape^''. Using a perturbative approach, it was found 
that even a small anisotropy stabilizes the exotic multi- 
quanta states which can therefore exist at lower Maki 
parameter (any > 1-8) than predicted by the ideal- 
ized isotropic models. According to this prediction such 



2 



states arc likely to oceur in any real anisotropic Pauli 
limited superconductor. More specifically in the tetrag- 
onal symmetry, 3 scenarios are possible for the FFLO 
state: a) Maximal FFLO modulation along the field with 
zero Landau level state, b) Highest Landau level modula- 
tion in the plane perpendicular to the field and no FFLO 
modulation, and c) Both Landau level and FFLO modu- 
lations. This three scenarios correspond to the tetragonal 
symmetry and was derived within a single Landau level 
approximation, which is valid at large field. It may thus 
be relevant to explain the observation of two high-field 
and low temperature phases of CeCoIns, which exhibit 
contrasted behaviors under distinct magnetic field ori- 
entations (inside or perpendicular to the Celna planes). 
Nevertheless the shape of the Hc{T) transition line is 
far from sufficient to establish a clear correspondence be- 
tween one phase and a particular class of solutions among 
the three a)-c) possibilities. It is thus necessary to gain 
complementary informations to determine which scenario 
among a)-c) is actually realized. As natural prccussors of 
the transition, the fluctuations in the normal state pro- 
vide informations about the superconducting state itself. 
We shall show here that fluctuations enable to detect the 
presence of a FFLO state in both 2D and 3D supercon- 
ductors, and allow to discriminate between the various 
a)-c) scenarios in the tetragonal 3D case. We concentrate 
on the region near the tricritical point (T*, H*) which is 
the meeting point of the three transition lines separat- 
ing respectively the normal state, the uniform and the 
modulated superconducting states^'"^^. 

In this paper, we evaluate the fluctuation induced 
magnetization near the FFLO transition in both 2D 
and 3D anisotropic superconductors using the modified 
Ginzburg-Landau (MGL) functional^^"^^. Previously we 
calculated the fiuctuational specific heat and conduc- 
tance near the pure FFLO transition in the absence of 
orbital cffect^^. Our motivation was to establish a re- 
lation between the topology of the lowest energy fiuctu- 
ation modes and the divergencies of the physical prop- 
erties at the FFLO transition. In the isotropic model, 
those divergencies arc very different than the standard 
BCS ones since the topologies of the degenerate FFLO 
and BCS modes differ fundamentally. Unfortunalely, in 
the anisotropic models, this degeneracy is lifted and the 
topologies of FFLO and BCS modes become quite simi- 
lar, thereby leading to less contrasted behaviors. 

In the two-dimensional case, we also show that the ra- 
tio between the paraconductivities along {axx) and per- 
pendicular {(Jyy) to an in-plane applied magnetic field 
H = Hbx is drastically enhanced near the FFLO tran- 
sition, in comparison to the one near the uniform BCS 
transition. Moreover we demonstate that the high-field 
fiuctuational magnetization of thin films may oscillate be- 
tween positive (paramagnetism) and negative (diamag- 
netism) values. These oscillations originate from the 
competition between orbital and paramagnetic effects 
which tend to promote respectively Landau level mod- 
ulation and FFLO modulation. Being precursors of the 



Meisncr or Abrikosov lattice state, the superconducting 
fluctuations are usually diamagnctic. Therefore the para- 
magnetism predicted here is a hallmark of the unconven- 
tional FFLO state. At lower field, the magnetization is 
suppressed near the FFLO transition in comparison to 
the BCS case. These features should also pertain in the 
case of layered 2D compounds like A-(BETS)2FeCl4 or 
K-(BEDT-TTF)2Cu(NCS)2. 

In 3D superconductors under high magnetic field, these 
oscillations are blurred out when scenario a) is real- 
ized whereas they pertain when scenario b) takes place, 
thereby providing an experimental test to distinghish 
among the various possible structures of the order pa- 
rameter described in Ref.^^. Experimentally, the super- 
conducting fluctuations in CeCoIns have been investi- 
gated far above Tc and under low fields^'^'^'*. Here we 
suggest similar measurements near under strong mag- 
netic field and near the FFLO critical temperature. 

The paper is organized as follows. In Sec II, we present 
the MGL formalism which includes higher order deriva- 
tive of the order parameter than the standard Ginzburg- 
Landau functional. Such an extension is necessary to 
handle the nonuniform FFLO state. In Sec III, we anal- 
yse the case of thin superconducting films under in-plane 
magnetic field and predict a strong dependence of the 
conductance upon the mutual orientation of the current 
fiow and magnetic field. We also derive the fluctuation 
magnetization induced by a tilted magnetic field point- 
ing out of the film plane. In Sec IV, we discuss the 
3D anisotropic compounds with emphasis on the fluc- 
tuation magnetization. In appendix, we provide a de- 
tained derivation of the Ginzburg-Levanyuk criterion for 
the FFLO transition in order to justify the Gaussian ap- 
proximation used in the whole paper. 



II. FORMALISM 

Here we introduce the modified Ginzburg-Landau 
(MGL) free-energy functional and detail how the fluctu- 
ation induced properties can be obtained from the spec- 
trum of the fluctuations. This approach is valid near the 
tricritical point. The location of the tricritical point in 
the temperature (T)-magnetic field {H) phase diagram 
depends on microscopic details of the model like the con- 
centration and the type of impurities, the crystal and the 
order parameter symmetries. Nevertheless the formula 
of this section are valid generically around the tricritical 
point, independently of its precise location. Note that 
for a clean s-wave superconductor the tricritical point is 
located at T* = OMTco,ij,bH* = 1.07fcsTco, T^o being 
the critical temperature in the absence of Zeeman effect 
and hb the Bohr magneton^'^^. 
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A. Modified Ginzburg-Landau functional 



Fluctuation free energy in the absence of orbital 
effect 



At the vicinity of the tricritical point {T*,H*), the 
FFLO transition can be described by a Ginzburg- 
Landau Uke approach since the order parameter ^'(r) 
and its spacial gradients are small^^. The correspond- 
ing free-energy functional, called hereafter the modi- 
fied Ginzburg-Landau (MGL) functional, diff'ers from 
the usual Ginzburg-Landau functional by the presence 
of higher order derivatives of the order parameter *(r). 
This is related to the fact that the FFLO phase cor- 
responds to a non-uniform groundstate. This MGL 
functional may be directly derived from the microscopic 
BCS theory for clean isotropic s-wave superconductors 
with Zeeman interaction'^'', and allows extention to con- 
ventional and unconventional singlet superconductors 
in the presence of paramagnetic, orbital and impurity 
effects20>2i. The quadratic terms (|^'(r)|^ |V*(r)| , 
etc..) of the MGL describe the free dynamics of the or- 
der parameter while higher-order terms (|^'(r)|^, |^'(r)|^ 
etc..) account for the interactions. In this paper, we 
analyse the fluctuations within the Gaussian approxima- 
tion whereby only the quadratic terms are retained^^. 
In the appendix, we shall justify this approximation by 
demonstrating that the interaction terms are small in the 
experimentally relevant range of parameter Tco/Ep <C 1, 
Ep being the Fermi energy. 

Within the Gaussian approximation, the MGL func- 
tional writes 



In the absence of magnetic field, the order parameter 
can be expanded in plane waves as \l/(r) = J2k'^^^^^''^ ■ 
In this Fourier representation, the free-energy Eq. (1) 
can be rewritten as Hq [^] = X^k^k |^k|^ where 
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Sk = a- Qiki +^ijkik 



(2) 



describes the spectrum of the decoupled fluctuation 
modes ^fk- The partition function Z = Tr(e~^'^/'^^"^) is 
obtained by tracing over all the possible values of these 
modes v^k which reduces to an infinite product of Gaus- 
sian integrals: 



z = n 

k 



(3) 



The corresponding thermodynamical free energy (per 
unit volume) F = —ksThiZ is given by^^ 



F = kBT 



J (2^ 



In 



Ek 



TvksT 



(4) 



In the isotropic case {gi = g and = 7), the spectrum 
Ek = a — -|- 7k^ can be exactly rewritten as 



Sk = T + (k^ - ql)"^ 



(5) 



where Qq = g/2j and t = a — g'^/4'j. The FFLO transi- 
tion arises at the critical temperature {H) defined by 
T = 0, namely at: 



T, (H) = % (H) 



47a' 



(6) 



Hg [*] = Jdr(a 1*1' - gi | A*|' + 7ij I A A*l') , 

a) 

where the summation over the index i = x, y, z is im- 
plied, and a = a (t — , Tc {H) being the critical tem- 
perature for the uniform superconducting second order 
phase transition. The components Ai of the vector po- 
tential enter the MGL through the covariant derivatives 
Di = di + 2ieAi/fi, while the coefficients a,gi and 7ij 
are functions of the temperature and field. The detailled 
microscopic expressions of these coefficients as functions 
of T and H can be found in^^~^^ for various crystal lat- 
tices and in presence of paramagnetic, orbital and impu- 
rity effects. As a salient and common feature of these 
functionals^'^"^'^, coefficients gi change sign at the tricrit- 
ical point, thereby inducing an inhomogeneous supercon- 
ducting phase when gi > 0, namely at low temperature 
and high field {H/T > H*/T*). Then the FFLO crit- 
ical temperature Tc (H) is larger than Tc {H) implying 
a transition between the normal state and a nonuniform 
superconducting state. Note that the idealized isotropic 
form of this functional corresponds to gi = g and jij = 7. 



which is higher than the critical temperature for the tran- 
sition towards a uniform superconductor Tc{H). In the 
normal state T > Tc (H), Eq. (5) makes apparent that 
the lowest energy fluctuation modes are degenerate and 
located around the sphere k^ = in reciprocal space. 



C. Fluctuation free energy with orbital effect 

Here we consider Eq. (1) in the isotropic case (gi = g 
and jij = 7) in presence of the orbital effect associ- 
ated with a magnetic field H = Hbz- Within the gauge 
Ax — and Ay = xH, the order parameter is conve- 
niently represented in terms of Landau wavefunctions. 
This follows from the fact that the mean-field equation 
SH/S'i!* = reads 



(a + .gD^ + 7D*)*(.T, y, z) = y, z) 



(7) 



Owing to the translational invariance along y and z, the 
momenta ky and kz are good quantum numbers. In the 

absence of the fourth-order derivative (7 = 0), the so- 
lutions of Eq. (7) are well known'*^, being the standard 
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Landau wavcfunctions ^'(r) = ^n,ky,k^ = /(a;) 6**^"^+*'^^ 
The equation for f{x) 



2e. 



Hx] f + ia-gki)f = Ef{x) (8) 



is similar to the harmonic oscillator equation with "in- 
verse mass" —g and the frequency lOc = —AeHg/H"^. 
Therefore the Landau levels for (7 = 0) are given by^^ 



En{kz) = a-gkl + fiwc(n + -) 



(9) 



where n = 0, 1, 2, .. while kz is a continuous wavevector. 

We now solve Eq. (7) in the presence of the fourth- 
order derivative (7 7^ 0) which is a hallmark of the 
MGL functional and FFLO state. We first show that 
the Landau wavefunctions ^n,ky,k^{''^) (eigenstates of 
a + .gD^) arc still eigenstates of the differential opera- 
tor a + gD^ + 7D^. Indeed introducing the quantized 
wavevector Qn (n = 0,l,2,...) defined by 



Qn = —y^+n) 



AeH 



one obtains immediately that 

5D^$n,fe„fe. = -g {kl + Q^) 

and therefore^ ^ 

with the fluctuation spectrum 



(10) 

(11) 

(12) 

(13) 
(14) 



Note that En{kz) is still degenerate, being independent 
of ky. Finally the free energy functional reduces to the 
sum 



(15) 



over decoupled modes $n.fe„.fe,- The partition function 
Z = Tr(e~^'^/'°^"^) is obtained by tracing over all the 
possible values of these modes ^n,ky,k^ and finally re- 
duces to an infinite product of Gaussian integrals: 



TrksT 



n,k},k^ En{kz) 



(16) 



The corresponding thcrmodynamical free energy (per 
unit volume) F = — ksT In Z is given by 



^ if , ^ f dk, ^ 
F = — keTy / — n 



Enjk,) 
TrksT ' 



(17) 



where the prcfactor H/^q accounts for the degeneracy 
of each Landau level at given n and kz^"^. Here $0 = 



h/2e is the superconducting flux quantum. In the two- 
dimensional case the quantum number kz is irrelevant 
and the average free energy per unit surface is given by 



F = — fcsTVln 
$0 ^ 

" n=0 



Enjkz = 0) 

irksT 



(18) 



The magnetization along the z-axis is simply given by 
M = -dF/dH. 



D. Transport 

Besides thermodynamics, the fluctuations also affect 
the transport properties. A standard procedure con- 
sists in using the time-dependent Ginzburg-Landau equa- 
tion to obtain the current-current correlator at different 
times^^'^^. Then a general formula for the paraconduc- 
tivity tensor 

AH J {2w)d 4 

can be obtained within the Kubo formalism^*. It should 

be noticed that the momentum dependence of the veloc- 
ity component Wki = de^/dki along the i— axis is very 
different from its usual form owing to the presence of the 
'Jijkfkj terms in the dispersion relation Eq.(2). In par- 
ticular, Wki are no longer linear combinaisons of the mo- 
mentum components ki. The formula Eq.(19) provides 
the so-called classical Aslamasov-Larkin contribution to 
the paraconductivity. It is well known that the quantum 
Maki-Thomson contribution can be important especially 
in the two-dimensional case^'*. In this paper, we study 
the FFLO transition at the vicinity of the tricritical point 
(T*, H*) where the strong pair-breaking mechanism sup- 
presses the Maki-Thomson contribution. 



III. SUPERCONDUCTING THIN FILMS 

Here we consider thin superconducting films where the 
FFLO state is realized due to the paramagnetic effect of 
an in- plane magnetic field H = H^^Bx- Wc assume that 
the lattice has the square symmetry and thus equiva- 
lent properties in the x and y directions in the absence 
of field. In Sec III. A, we first neglect the orbital effect 
associated with H = e^; and calculate the paraconduc- 
tivity o"xx = cTj/y Then we treat perturbatively the or- 
bital effect associated with H. = H^^Bx, and find that the 
paraconductivity axx measured along the applied mag- 
netic field differs from the one {(Jyy) measured perpen- 
dicular to the field (Sec. III.B). Finally, we also discuss 
the effect of a perpendicular magnetic field H = Hbz 
which induces magnetization oscillations between dia- 
magnetism and paramagnetism (Sec. III.C). This behav- 
ior is in sharp contrast with the usual fluctuation induced 
diamagnetism predicted and observed close to the BCS 
transition^^'^^. 
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FIG. 1: Minima of the fluctuation spectrum for the two- 
dimensional model with square symmetry. The filled (open) 
ellipses indicate the locations of the lowest energy fluctuations 

for ?7 > 1 (77 < 1). Note that in the isotropic model {rj = 1), 
there is an infinity of degenerate minima located on the circle 
(solid line). 



A. Pure pEiramagnetic limit 

In the pure paramagnetic limit, the in-plane magnetic 
field H = H^^Bx only acts on the spins. Hence one can set 
= in the MGL functional Eq.(l) and use the spec- 
trum Eq.(2) to describe the superconducting fluctuations 
in the normal state. By increasing iJ|| and lowering the 
temperature, one may tune the thin film near the tricrit- 
ical point where g = 0. The parameter g may change 
sign, thereby indicating an instability towards the inho- 
mogeneous FFLO state when g is positive. 

Furthermore the square symmetry imphes (^a; = gy = g, 
Ixx = lvv=l and ^xy = lyx = Eq.(2). Isotropy is 

only restored when r] = 1 which corresponds to a circular 
Fermi surface^^. Then the low energy modes are located 
around the whole circle defined by = = 5/27 in 
reciprocal space, according to Eq.(5). In the general case 
(rj 1), the quartic terms Jijkfk'^ introduce a non triv- 
ial anisotropy, and the lowest energy are realized around 
isolated points of the reciprocal state (Fig.l). 

The paraconductivity (t^x = o'yy near the FFLO tran- 
sition (g > 0) is obtained from Eq.(19) where the integral 
extends over the whole reciprocal space. Nevertheless 
due to the denominator in Eq.(19), the main contri- 
bution comes from the lowest energy fluctuation modes. 

We shall flrst identify these isolated minima kmin of the 
energy by solving 5ek/9fca,(kmin) = 5£k/9/cy(kmin) = 0. 
The location of these minima differs in the cases r] > 
1 and rj < 1 respectively (Fig.l). In the case rj > I, 
the spectrum £k has four minima located along the x 
and y axis, namely at points A± {±qo, 0) and B± (0, ±qo) 
of the reciprocal space (Fig.l). Near the minima A± 
located along the x axis, we find that the spectrum can 



be approximated by 

slf^^ ^T + 2g{kxTqof + g{v-'^)kl (20) 



with T = a — = a (T — Tg). The expression for the 

FFLO critical temperature Tc coincides with Eq.(6) as 
long as > 1. 

Similarly, the spectrum around the minima 
B± (0, ±qo) located along the y axis is given by 



ei^^^ ^T + g{r,-l)kl + 2g{kyTqof. 



(21) 



The next step consists in evaluating the generalized 
velocities around each minima. For instance around 
we obtain 



de 



dkx 



= 4.g{kx - qo) ■ 



=2g{r,-l)ky. 



dk 



(22) 
(23) 



■V 



Finally, we evaluate the integral Eq.(19) around the min- 
imum j4_|-. 

Reproducing this steps for the other minima A- , 
and B- and summing the contributions from the four 
minima, one obtains the paraconductivity 



^xx — ^ 



e^akB 1 



vv 



4V2H \T-T, 



(24) 



which diverges at the FFLO transition. 

The second case r] < 1 can be treated along the 
same line of reasoning, albeit the four degenerate; minima 



k' = (k' 

mill V n 
with (fc^ 



mm X ' mm y 



) are now located on the diagonals 
(kLinyf = 90/2 (Fig.l). The spectrum 



can be expanded as 

Ek ~ T — Xij {ki — fcjjjjjjj) (fcj — fcjjjjjjj-) , (25) 

i,j—x,y 

around any of those four minima denoted A'^ , A'_ , B'^ 
and B'_ (Fig.l), with r = a-g^j^-f (1 -f 77) = a (T - Tc). 
For instance, we find 



2fl 
1 + 



and \xy — ^yx 



2m_ 

1 + r? 



(26) 



aroimd ^'^((jo/\/2, (/o/v^). Diagonalization of the ten- 
sor Ay leads to the eigenvalues \x = 2g and Ay = 
2g (1 — 77) / (1 -|- 7]) along the principal axis X and Y. 
Finally summation over the four minima (Fig.l) yields: 



'vv 



1 + V 



2fi {l + r]y\l-v\T-Tc 



(27) 



Note that in this regime 77 < 1, the expression for the 
FFLO critical temperature Tc differs slightly from Eq.(6). 

The above expressions Eqs. (24,27) both diverge for 
77 ^ 1, which indicates stronger Gaussian fluctuations 
in the isotropic model^^. 
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B. Orbital effect associated with an in-plane 
magnetic field 

Wc now take into account the orbital effect associated 
witli the in- plane magnetic field (H = H\]^ex)^ and show 
that it breaks the square symmetry, inducing distinct 
paraconductivities along x and y directions. 

In the case of thin films with a strong confinement 
in the z direction, this orbital effect is small. Then it is 
still possible to describe the fluctuations by the spectrum 
Eq.(2) with iJii -dependent coefficients^^, namely by 



: a 



(Ml 

12<1>2 



+ 7 



+ 



80$4 



+ 27(r?-l)fc^fc2+ M_(^^fc2^3^fc^), (28) 

where k"^ = + ky, d is the width of the film along the 
^-axis, and $o = h/^e is the superconducting quantum 
of flux. Owing to the smallness of the dimensionless pa- 
rameter (ffy^d/go'I'o)^, it is still possible to use Eq.(19) 
in order to evaluate the paraconductivity tensor follow- 
ing the same procedure than in the previous section III A. 
Here ^ is the superconducting coherence length. 

Due to the field (iJ|| ) dependence of the coefficients in 
Eq.(28), the minima of £k are displaced (in comparison 
to the case H^^^d/qo^o = shown in Fig. 1) according 
to {r] > 1): 



(29) 
(30) 




Moreover the square symmetry is broken since the critical 
temperature associated with modulation A± differs from 
the one for B±. It happens that the FFLO modulation 
occurs along the field (points A±) for ry < 3. In contrast 
for > 3 the modulation occurs along the y-axis (points 
B±) which is perpendicular to the applied field^^. 

We concentrate on the case for 1 < r/ < 3 wherein 
the order parameter is modulated along the field H = 
-ff||ej.. The paraconductivity comes from the contribu- 
tions around the points and A_ . Then the paracon- 
ductivity CTxx measured along the field differs from the 
one ayy measured in the perpendicular direction. As a 
main result of this section, the ratio (Jxx/'^yy contains 
a contribution which diverges at the FFLO transition 
(ff>0) 



1 + :?^ 



H\\d 



6g \ $0 



T] — 1 2jT T] +1 



(31) 

This term produces a strong enhancement of o^xl'^yy 
when the transition is approached, namely when t = 



a{T — Tc) 0+. This is in sharp contrast with the 
transition towards a uniform BCS state. There the ratio 
of the paraconductivities 



cr, 



= 1 - 



7(r?-3) (li\\d 



VV / BCS 



eg 



(32) 



does not contain any divergent term at the BCS transi- 
tion {g < 0). Such an enhancement of <7xx/<^yy may serve 
as an experimental signature of the FFLO state. This 
property is reminiscent of the recent finding that critical 
current oscillates as a function of magnetic orientation in 
anisotropic 2D films^^. 



C. Orbital effect in a perpendicular magnetic field 

Finally, we discuss the effect of a perpendicular compo- 
nent H±ez when the field H = H\\ex + H±ez is tilted out 
of the film plane. The perpendicular component H±ez 
quantizes the in-plane motion of the fluctuating Cooper 
pairs, and induces a finite magnetization. This effect is 
larger than the orbital motion associated with the in- 
plane part of the field. We therefore neglect the later 
(which is correct for H\\£^d <^ <I>o) and use the MGL fimc- 
tional Eq.(l) within the gauge Ax = and Ay = xH±, 
like in section II. C. 

In the anisotropic case, the eigenmodes of this func- 
tional are not known exactly for finite H± precluding 
an analytical evaluation of the magnetization. Never- 
theless the isotropic model already exhibits oscillations 
of the fiuctuation magnetization. Recently the fluctua- 
tional magnetization (persistent current) of small rings 
made of a FFLO superconductor was obtained within 
the framework of the isotropic model introduced in Sec. 
II. C. In the following, wc derive a simple formula for the 
magnetization in the simpler planar geometry of super- 
conducting films^°. 

The fluctuation spectrum is given by Eq.(14) and the 
free energy per unit surface by Eq. (18) with H = H±. 
The particular form Eq.(14) of the spectrum enables both 
degeneracies between the Landau levels (£„ = En+i) 
and commensurability effects between the wavevectors 
Qn and go- 

Single mode (high fields ) approximation: For large per- 
pendicular field, namely H±/^o ^ v^T/7, the Landau 
levels are well separated from each others and the main 
contribution to the free energy Eq.(18) comes either from 
the single level with minimal energy En , or from two lev- 
els when a degeneracy (En = -Eji+i) occurs. 

Let us first consider the non-degenerate case. Then the 
free energy is simply given by the single level contribution 



$0 T^ksT 



(33) 



and the corresponding orbital magnetization (per unit 
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surface) 



(8n + 4)fcBT/ -f{Ql~q^o) 



(34) 



is highly nonUnear since the prefactor of H± depends 
strongly on the field and on temperature. Importantly 
the magnetization may change sign due to the presence 
of the factor — Qq in the numerator. In order to make 
more transparent the formula Eq.(34), one may intro- 
duce the field-dependent temperature Tcn{H) where the 
denominator vanishes: 



a{T-T,n) = T + 'r{Qi-q^ 



,2\2 



(35) 



This relation defines the second-order transition line 
Tc"\h) between the normal and the modulated super- 
conducting state described by the n-th Landau level. We 
also define the points like A,C,E (Fig. 2) along this tran- 
sition line where the numerator vanishes since Qn = Qo- 
Those points are also located on the second-order transi- 
tion line Tcp{H) between the normal and the FFLO su- 
perconducting state calculated in the pure paramagnetic 
limit. In the normal state, the orbital magnetization can 
be therefore reexpressed as 



Mr,. = - 



{8n + 4)kBT'y fQ 



Qo 



(36) 



the tricritical point is approached, i.e. when g — > 
(Fig.2)22. On the FFLO side (region g > in Fig.2) and 
for a given magnetic field, the system becomes a FFLO 
superconductor before the divergency develops because 

Tc > T,. 

Conclusion: At high fields H±, the magnetization near 
the FFLO transition line oscillates between sizeable dia- 
magnetism and paramagnetism as the transitions be- 
tween successive Landau levels are realized. In the low 
field limit, these transitions become very close and the 
oscillations average themselves leading to a cancellation 
of the linear response and a suppression of the diver- 
gency. This situation is in strong contrast with the stan- 
dard BCS case where the magnetization is always dia- 
magnetic. We suggest to perform magnetization mea- 
surements in thin films near an expected FFLO transi- 
tion. The suppression of the fluctuational magnetization 
at low perpendicular field and the restoration of sizeable 
oscillations between para- and diamagnetism at higher 
perpendicular field should be a strong indication for the 
FFLO state in quasi-two dimensional compounds. 



This 2D magnetization is diamagnetic when > qQ 
and paramagnetic when Q"^ < q^ (Fig.2). In contrast, 
the fluctuation magnetization is always diamagnetic in 
the BCS case. However M„ follows a similar power law 
(T - ri"))-i and is on the same order of magnitude 
than the BCS magnetization^^. Consequently we expect 
that the oscillations between diamagnetism and param- 
agnetism should be measurable in thin films of FFLO su- 
perconductors. This single mode approximation breaks 
down when the n-th and {n + l)-th Landau levels are de- 
generate, namely when En = En+i- Then the two levels 
must be included together in the free energy, whereas the 
other Landau levels are still far in energy and can be ne- 
glected safely. The resulting magnetization M„ -I- M„+i 
is slightly diamagnetic at degeneracy. 

Continuum (low fields) approximation: The single 
mode approxi matio n breaks down in the weak field limit 
{Hj_/^Q <^ \J t/^) where the Landau level separation 
becomes so small that all the levels have to be taken into 
account. This situation corresponds to a magnetic field 
which is slightly tilted out of the film plane. In the case 
of a uniform BCS superconductor, the standard result for 
the magnetization is^^: 



3 a{T-T,) 



(37) 



which is diamagnetic {g < 0) and diverges at the BCS 

critical temperature T^. (H) for the second order phase 
transition. This diamagnetic response is suppressed when 



IV. ANISOTROPIC 3D SUPERCONDUCTORS 



It is commonly believed that the FFLO state in 
CeCoIus corresponds to a modulation along the applied 
magnetic field. Nevertheless it was argued recently that 
this situation is unlikely to happen for arbitrary field ori- 
entations when the tetragonal anisotropy of CeCoIus is 
properly taken into account. Apparently if the order pa- 
rameter modulation is along the field for H ± c (resp. 
H II c), then the modulation is likely to be perpendicu- 
lar to the field for H || c (resp. H ± c)^^. Here we in- 
vestigate the FFLO fluctuations in anisotropic 3D com- 
pounds, building upon the various mean- field scenarios 
reported in Ref'^^. We evaluate the fluctuational mag- 
netization M along the magnetic field H = Hez- In 
particular, we demonstrate below that the magnetiza- 
tion oscillates between sizeable diamagnetism and para- 
magnetism when the modulation is perpendicular to the 
field (Landau level like). Those oscillations are the 3D 
counterparts of the ones predicted in the previous section 
for superconducting films. In contrast the magnetization 
is shown to be strongly suppressed when the modulation 
occurs along the field (FFLO like modulation). In the 3D 
case, magnetization measurements therefore provide an 
experimental tool to discriminate between the two possi- 
ble order parameter structures uncovered in Ref^^. 
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A. Mean Held 

We start by a short reviewing of the mean-field prop- 
erties of the functional 



if = a|*| 



E 

i—x,y,z 



<?|A*|'+7 



i—x,y.z 



\Dtm 



(38) 



consistent with the tetragonal symmetry of CeCoIn5. 
The terms e^, Sx and e describe nontrivial (namely dif- 
ferent from a simple elliptical) anisotropy-'^^. Note that 
the cubic symmetry corresponds to e^; = e and = 0. 
It was shown that two kinds of modulated superconduct- 
ing states (scenarios a) and b) mentioned above in the 
general introduction) are the most likely to occur when 
anisotropics arc properly taken into account. 

The class a) of solutions corresponds to order parame- 
ters modulated along the field with characteristic FFLO 
wave-vector and in the n = Landau level. Following 
the mean field analysis of Ref^^, we write the fluctuation 
spectrum as 



2eH 



(39) 



which indicates an instability towards finite modulation 

along the z axis (magnetic field). This is similar than 
Eq.(14) except that the Landau index is fixed n = and 
7a is a renormalized parameter (specific to this class a) 
of solutions) which depends on e^, Sx and e. 

In the class b) of solutions, the modulation occurs in 
the plane perpendicular to the field and is described by 
a higher {n > 0) Landau level. Following the mean field 
analysis of Rcf^^, we write the fluctuation spectrum as 



En{k,) =T + jbiQi - Qo) + 9bk 



(40) 



This spectrum differs from Eq. (39) since the kinetic en- 
ergy term gtfc^ favors kz — 0. Furthermore a finite Lan- 
dau index is possible, since the energy is also minimized 
by choosing n such as — 9o 0- brief, the lowest 
energy fiuctuations in the normal state ressemble the su- 
perconducting groundstate which is modulated perpen- 
dicularly to the field. Note that here 75 and gi, are also 
renormalized parameters which are specific to the class b) 
of solutions and depend on e^, Sx and e in a complicated 



manner 
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B. Fluctuation magnetization 

Here we evaluate the magnetization induced by the 
FFLO fiuctuations taking into account the intrinsic 



anisotropy present in 3D compounds. The FFLO transi- 
tion might happen under low or strong field, depending 
on the underlaying microscopic mechanism. For instance, 
in the rare earth magnetic superconductor ErRh4B4 a 
small field is sufficient to polarize the internal moments, 
and the FFLO transition is thus expected at low applied 
magnetic field^. Here we treat the case of the FFLO 
transition occuring under strong magnetic field which is 
relevant for the case of the heavy fermion superconduc- 
tor CeCoIus. Using a single Landau level approximation, 
we demonstrate that the magnetization exhibits qualita- 
tively distinct behaviors depending on the class of solu- 
tions. 

FFLO-like modulation along the field, characterized by 
a finite wave-vector qo and Landau index n = ( scenario 
a) discussed in the introduction). The 3D density of free 
energy is given by the integral 



H 



knT— 



dkz . E„^o{kz) 
■ m 



$0 J 27r 



(41) 



where the energy En=o{kz) is given by Eq.(39). Hence 
the most divergent part of the orbital magnetization (per 
unit volume) M = —dF/dH is given by 



M 



2kBT 



{/.: 



dkz 



ja{2eH/H+kl 



T + 7a (2ei//fi+fc2-Q2)^ 



H, (42) 



where $0 = h/{2e). Since the numerator of the inte- 
grand cancels and changes sign as a function of fc^, one 
expects a strong suppression of the fluctuation magneti- 
zation compared to the uniform BCS case wherein such 
a cancellation does not occur. Indeed the magnetiza- 
tion M diverges logarithmically at the FFLO transition 
which is less divergent than the t"^/^ law predicted in the 
standard BCS case. Therefore the presence of a genuine 
FFLO state should be detected as a suppression of the 
fluctuation diamagnetism observed near the BCS tran- 
sition. In comparison with the 2D case, the oscillations 
between paramagnetism and diamagnetism predicted in 
the previous section arc blurred out by the dispersion 
over the momentum kz along the fleld. 

Landau level modulation perpendicular to the field (sce- 
nario h) discussed in the introduction). The dispersion 
of the fluctuations Eq.(40) now favors the absence of 
modulation along the z axis in contrast to the spectrum 
Eq.(39). Upon increasing the parameter q^, the lowest 
energy Landau level is successively n = 0, then n = 1 
etc... Near the BCS transition {g > 0), the fluctua- 
tions induce diamagnetism and a lowering of the critical 
field Hc2{T) below the purely paramagnetic critical field 
Hp{T) at the same temperature (Fig. 2). When the n-th 
Landau level is realized and when all the other Landau 
levels are distant in energy, one can single out the con- 
tribution of this main level to the density of free energy 



°° dK, T + 76(Q^ - qlf + gbkl ^^^^ 



27r 



TrksT 
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Writing the orbital magnetization as 



$2 



dkz 



IbiQl - ql) 



H±, (44) 



shows that the 2D oscillations arc no longer suppressed 
by the integration over since the numerator is inde- 
pendent of kz- Calculating the integral shows that the 
magnetization diverges as 



M„ = - 



{4n + 2)kBT 



lb 



Ql 



ql 



$2 



[agbY^'^ (T - Ten) 



1/2 



(45) 



with the same power law than in the BCS transition of 
3D superconductors^^. Unlike the BCS case, this fluctua- 
tion magnetization changes sign being diamagnetic when 
<5n > 9o (arcs BC, DE in Fig. 2) and paramagnetic when 
Ql < ql (arcs AB, CD). 

In brief, the magnetization is sizeable and oscillates be- 
tween para- and diamagnetism when the superconduct- 
ing order parameter is modulated perpendicularly to the 
field, whereas it is strongly suppressed when the order 
parameter is modulated along the field. Therefore mag- 
netization measurements may serve as a test to discrim- 
inate between FFLO and Landau level like modulations 
in 3D anisotropic superconductors. 



V. CONCLUSION. 

We investigated the conductivity and the orbital mag- 
netization associated with superconducting fluctuations 
above the FFLO critical temperature or field. Both in 
2D and 3D models, we shown that these properties dif- 
fer considerably than their counterparts at the vicinity 
of a standard BCS transition towards an homogeneous 
superconducting state, thereby providing an experimen- 
tal tool to detect the inhomogeneous state. First, the 
paraconductivity of thin superconducting films exhibits 
a strong anisotropy when measured parallel or perpen- 
dicular to the FFLO modulation. Second, the orbital 
magnetization oscillates between diamagnetic and para- 
magnetic behaviors at high fields, and is strongly sup- 
pressed at low fields, whereas the uniform BCS state 
always induces diamagnetic fluctuations above Tc- We 
suggest performing magnetization and conductance mea- 
surements along the FFLO transition line in compounds 
where the FFLO state has been recently reported. In 
2D organic superconductors^' the magnetization os- 
cillations should be even more pronounced than in the 
3D magnetic superconductors (ErRh4B4 , see"^) or in the 
case of the anisotropic 3D heavy fermion superconduc- 
tor CeCoIns*'"^. It was recently shown that CeCoIns 
has quasi-2D Fermi-surface sheets coexisting with a 3D 




FIG. 2: Schematic ficld-tcmporature {H, T) phase diagram 
showing the cascade of Landau levels. The thick dashed curve 
represents the critical field in the absence of orbital effect. 
In presence of orbital effect, the critical field is reduced and 
described by the thin dashed curves which corresponds re- 
spectively to the n = 0, n = 1 and n = 2 Landau levels. The 
solid line represents schematically the expected transition hne 
between the normal and the superconducting states. The fluc- 
tuations are diamagnetic between zero field {H = 0, Tc) and 
B, paramagnetic near the arc AB, then again diamagnetic 
near the arc BC, etc. ..The Landau levels are degenerate at 
points B,D, etc... This schematic picture is relevant for both 
2D FFLO superconductors and for 3D ones where the Lan- 
dau level modulation is realized (scenario b) evoked in the 
introduction). 



anisotropic Fermi surface^**. Nevertheless due to strong 
hybridization, superconductivity in CeCoIns is likely to 
be described by a single order parameter. The com- 
plex structure of the Fermi surface should only modify 
the expressions of the coefficients in the MGL functional 
as functions of the microscopic parameters of CeCoIns. 
Moreover, the fact that the Lande factor is anisotropic 
leads to different Maki parameters depending on the field 
orientation and may also shift the position of the tri- 
critical point®. Nevertheless the magnetization oscilla- 
tions/suppression predicted here are generic of the pres- 
ence of FFLO phase, and should pertain independently 
of the microscopic characteristics of CeCoIns. 

Finally, in the 3D case, we find that the absence of such 

oscillations reveals a FFLO state modulated along the 
field whereas presence of oscillations should be associated 
with a multiquanta Landau modulation perpendicular to 
the field. 
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VI. APPENDIX 



In this appendix, we address the validity of the Gaus- 
sian approximation used in this paper. When the tem- 
perature is sufficiently close to the critical one, interac- 
tions between the fluctuation modes become so strong 
that the Gaussian approximation breaks down. In order 
to quantify the range of temperature where this break- 
down occurs, we derive the Ginzburg-Levanyuk criterion 
for the FFLO transition^s.ss.ae. 



A. Ginzburg-Levanyuk criterion for the FFLO 
transition 

The full isotropic MGL functional H [*] = Hq [*] + 
Hint [^'] contains a quadratic part 



Hr 



• m = NM j 



dr 



2 , M |«2,T,|2' 



|5*r + r ^ * 



(46) 
(47) 



same point (H, T) of the phase diagram (the tricritical 
point) For this reason and since we are solely inter- 
ested in orders of magnitude here, we have denoted the 
coefficient of the l*]'^ term by the same g used for the 

For examples of phase transitions, one usually eval- 
uates only the |^'|^ terms^^. Here the situation of the 
FFLO transition is rather specific since on the line g = 
g = of the (T, H) phase diagram the coefficient of the 
l^"!^ term vanishes. Therefore one should evaluate the 
next interacting term, l^*! , for the regions near this line 
g = 0. Sufficiently far away from this lino g = (see the 
quantitative criterion below), one may simply evaluate 
the l^*! term. 

Using Wick theorem to evaluate Hi„t [^] = -ff^t [^] + 
hI^I [*], we find 



2gNd{0) 



E 

k,k' 



(49) 



for the \^\ term, and 



^^L1[*] = iV.(0)Ef^E(l*-l^)o)'K .0 

k \ k' / 

(50) 

for the l^*! term. Note that in this problem the form of 
the free field correlator 



r + ^^ik^-q^) 



2 • 



(51) 



is rather special due to the proximity of the FFLO tran- 
sition. 



which describes the free dynamics of the order parame- 
ter ^f, and non quadratic terms ifj„t[^'] which describe 
interactions^^^^ . All the results of this paper are derived 
within the Gaussian approximation which consists in us- 
ing Hg [^] as the free energy functional thereby neglect- 
ing completely Hint [^]- In the spirit of the original^^'^^ 
and textbook^"'' Ginzburg-Levanyuk criterion, we evalu- 
ate the interaction terms 



Hint m=Na{0) j dr 



(48) 



in order to compare them with Hq [^]- 

Wc have introduced the dimcnsionlcss coefficients a = 
{T-fa)/fe, T = {T-Tc)/Tc, and 5 to make apparent the 
order of magnitude of ech term in the MGL functional as 
a function of the the energy scales Tc and Ep- In partic- 
ular the dimensionless parameter g is of order one. We 
have also introduced the rf-dimensional electronic den- 
sity of states Nd{0) and the superconducting coherence 
length ^ = vp/Tc (we set H = 1). Here = in 
analogy with the transformation performed in Eq. (5). 
It is a rather particular property of the MGL functional 
that the coefficients of the and vanish at the 



B. Isotropic model 

Far from the tricritical point, namely when Ij ^ 
{Tc/EFf^'^~'^^/^^~'^\ the leading correction to the Gaus- 
sian behavior originates from the interaction term. 
Morover the fluctuations propagate with a quadratic dis- 
persion, and the correlator Eq.(51) can be approximated 
as 



^ 7TkBT,/NdiO) 
'0 r + 4qU^{k-qo) 
when evaluating the sum 
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El* 



rf/2~-l/2 



T.Ndm^ 



(52) 



(53) 



Using TcNdiO)^'^ ^ {Ef/Tc)'^~^, we find that the j^-]^ in- 
teraction terms are negligible in comparison to the Gaus- 
sian ones when the condition (Ginzburg-Levanyuk crite- 
rion) 



:d/3 / y= 



Ep 



2((i-l)/3 



(54) 
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is fuUfilled. The critical region width is larger than in 
the standard BCS case (r > {TjEpf for d = 3 and 
r ^ Tc/Ep for d=2) but it remains extremelly thin. 

Near the tricritical point, when g < 
{Tc/Epf^'^-'^y^^-'^\ the 1*1^ interaction becomes 
stronger than the l^*!^ one since this latter contribution 
is suppressed by the extremelly small prefactor g. In 
particular, along the g = g = line in the {H,T) 
diagram, the |V'I'|^ and \'^\* terms are totally absent 
from the functional^''. Therefore one should compute the 



mean value ( l^* 



with a purely quartic momentum 



dependence. Since 



r2 Nd{0)L'i 



~(d-4)/4 



(55) 



the condition to neglect the l^*] 
fluctuation modes is thus 



interaction between the 



T > 



y,^ . 4(d-l)/(6-d) 

E^J 



(56) 



The critical fluctuations are present in a larger region 
of the phase diagram than for BCS superconductivity^^. 
During the completion of this work, we became aware 
of Ref.'^*' where the Ginzburg-Levanyuk criterion is de- 
rived by evaluting exclusively the |\['| interaction term. 
We therefore obtain the same Ginzburg-Levanyuk cri- 
terion as in Ref.^*^ for the large g regime whereas our 
criteria differ when approaching the tricritical point. In 
spite of this discrepancy, both procedures lead to the 
same practical conclusion that the critical region re- 
mains extremelly thin and inaccessible for experimen- 
tal observations, because of the smallness of the ratio 
Tc/Ep ~ (10-2 - 10-3). 



Moreover the fluctuations propagate with a quadratic 
dispersion, and the correlator Eq.(51) can be approxi- 
mated as 



Evaluating the sum 



^ TTkBTc/NdjO) 

T + gk'^(,'^ 



(d-2)/2 



(57) 



(58) 



and using TaNd{0)^'^ = (Ep/Tc)'^-^, we find that the 
interaction terms are negligible in comparison to the 
Gaussian ones when the condition (Ginzburg-Levanyuk 
criterion) 



i2-d)/{4-d) 



Tc_ 



2(d-l)/(4-d) 



(59) 



is fullfilled. This Ginzburg-Levanyuk criterion is similar 
(same power of Tc/Ep) than the one encountered in the 
standard BCS case. 

Near the tricritical point, when g < 
{Tc/Epf'-'^-^'^/'-^-'^^ the lowest energy fluctuations 
are located around the origin of the reciprocal space 
and have a quartic dispersion like in the isotropic model 
studied above. The Ginzburg-Levanyuk criterion is thus 
again 



T > 

near the tricritical point 



Tc_ 



4(d-l)/(6- 



(60) 



C. Anisotropic model. 

We now derive the Ginzburg-Levanyuk criterion in the 
case of anistropic FFLO superconductors. The large g 
regime is modified in comparison to the isotropic case, 
since there the low energy fluctuations are located around 
few isolated points instead being spread over a large shell 
of radius go • 

Far from the tricritical point, namely when g ^ 



(Tc/Ef 



^2(d-l)/(6-<i) 



, the leading correction to the Gaus- 



sian behavior originates from the interaction term. 



D. Conclusion 

We have obtained that the size of the critical re- 
gion in FFLO superconductors is more extended than 
in the usual uniform superconductor case. Neverthe- 
less, for superconducting compounds, the ratio Tc/Ep ^ 
(10^2 — IQ-^) is small and the critical region thus re- 
mains hardly accessible for experimental observations, 
thereby supporting the Gaussian analysis performed in 
this paper. This fact makes very difficult the observa- 
tion of the phenomena (first-order transition) predicted 
by renormalization group studies'^^'^^. 



* Also at Institut Universitaire de France; Electronic ad- 
dress: a . bouzdineQcpmoh . u-bordeauxl . f r 

1 P. Fulde, and R. A. Ferrell. Phys. Rev. 135, A550 (1964). 

^ A. Larkin, and Y. Ovehinnikov. Sov. Phys. JETP 20, 762 
(1965). 



3 L.N. Bulaevskii et al, Adv. Phys. 34, 175 (1985). 
* R. Casalbuoni, and G. Nardulh. Rev. Mod. Phys. 76, 263 
(2004). 

^ Y. Mastuda, and H. Shimahara J. Phys. See. Jap. 76, 
051005 (2007). 



12 



® H. Radovan et al. Nature 425, 51 (2003). 

^ A. Bianchi et al. Phys. Rev. Lett. 89, 137002 (2002). 

* A. Bianchi et al. Phys. Rev. Lett. 91, 187004 (2003). 

^ S. Uji et al. Phys. Rev. Lett. 96, 157001 (2006). 
1" R. Lortz et al. Phys. Rev. Lett. 99, 187002 (2007). 
" L.W. Gruenberg and L. Gunther, Phys. Rev. Lett. 16, 996 
(1966). 

A.l. Buzdin and J. P. Brison, Phys. Lett. A 218, 359 (1996). 
" A.I. Buzdin and J.P. Brison, Europhys. Lett. 35, 707 

(1996) . 

" J.P. Brison et al, Physica C 250, 128 (1995). 

K. Yang, and S. Sondhi. Phys. Rev. B 57, 8566 (1998). 
^® A. Vorontsov, J. Sauls, and M. Graf. Phys. Rev. B 72, 

184501 (2005). 

^"^ D. Denisov, A. Buzdin, and H. Shimahara, Phys. Rev. B 
79, 064506 (2009). 

D. Saint- James, D. Sarma, and E.J. Thomas, Type II su- 
perconductivity (Pergamon Press, New- York, 1969). 
1^ A. I. Buzdin, and H. Kachkachi. Phys. Lett. A 225, 341 

(1997) . 

D. Agterberg, and K. Yang. J. Phys.: Condons. Matter 
13, 9259 (2001). 

M. Houzet and V. P. Mineev, Phys. Rev. B 74, 144522 
(2006). 

F. Konschelle, J. Cayssol, and A. I. Buzdin, Europhys. 
Lett. 79, 67001 (2007). 

R. Bel et al, Phys. Rev. Lett. 92, 217002 (2004). 

Y. Onose, Lu Li, C. Petrovic, N. P. Ong, Europhys. Lett. 

79, 17006 (2007). 

L. D. Landau, and E. M. Lifshitz. Course of theoret- 
ical physics, Volume 5 : Statistical Physics, 3rd edition 



(Butterworth-Heinemarm, 1980). 
^® M. Tinkharn, Introduction to Superconductivity (Mc- 
Graw Hill, second edition, 1996). 

W.J. Skocpol, and M. Tinkham, Rep. on Prog. Phys. 38, 
1049 (1975). 

^® A. Larkin, and A. Varlamov. The physics of super- 
conductors, Vol. I : Conventional and high-Tc supercon- 
ductors, chapter 3, Fluctuation phenomena in supercon- 
ductors (K.H. Bcnncmann and J.B. Ketterson, Springer- 
Verlag, 2003). 

A.I. Buzdin, Y. Matsuda, and T. Shibauchi. Europhys. 
Lett. 80, 67004 (2007). 
^° A. A. Zyuzin, A.Yu. Zyuzin, arXiv:0901.1562; A. A. Zyuzin, 
A.Yu. Zyuzin, JETP Lett. 88, 147 (2008). 
This useful property relies on the isotropy of the model 
treated in Sec. II. In the anisotropic case (Sec. Ill) 7D* is 
replaced by a more complicated operator 7ijDiD| whose 
eigenmodes differ from the ones of and are not known 
exactly. 

L. D. Landau, and E. M. Lifshitz. Course of theoreti- 
cal physics, Volume 3 : Quantum Mechanics, 3rd edition 
(Butterworth-Heinemarm, 2001). 
A. Schmid. Phys. Rev. 180, 527 (1969). 
A. Koitzsch et al., Phys. Rev. B 79, 075104 (2009). 
A. Lcvanyuk. Sov. Phys. JETP 36, 571 (1959). 

^® V. L. Ginzburg. Sov. Solid State 2, 61 (1960). 
S. A. Brazovskii. Sov. Phys. JETP 41, 85 (1975). 

^® D. Dahdovich, and K. Yang. Phys. Rev. Let. 93, 247002 
(2004). 




1 



